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Abstract: We consider one-dimensional translationally invariant quantum spin 
(or fermionic) lattices and prove a Mazur-type inequality bounding the time- 
averaged thermodynamic limit of a finite-temperature expectation of a spatio- 
temporal autocorrelation function of a local observable in terms of quasi-local 
conservation laws with open boundary conditions. Namely, the commutator be- 
tween the Hamiltonian and the conservation law of a finite chain may result in 
boundary terms only. No reference to techniques used in Suzuki's proof of Mazur 
bound is made (which strictly applies only to finite-size systems with exact con- 
servation laws), but Lieb-Robinson bounds and exponential clustering theorems 
of quasi-local C* quantum spin algebras are invoked instead. Our result has an 
important application in the transport theory of quantum spin chains, in partic- 
ular it provides rigorous non-trivial examples of positive finite-temperature spin 
Drude weight in the anisotropic Heisenberg XX Z spin 1/2 chain [Phys. Rev. 
Lett. 106, 217206 (2011)]. 

1. Introduction 

1.1. The problem. In 1969, Peter Mazur suggested |13j that the time-average 

A = lim(_i.oo /g lit' A{t') of a bounded observable A{t) can be bounded from 
below by means of exact conservation laws, i.e. observables Qk, k = 1,2..., 
satisfying (d/di)(5fc — and being mutually in involution, namely 

(^^),= lnjij/dt'(A(0)A(O>,>E^- (1) 

(•)^ is a thermal average at inverse temperature /?, and observables Qk have to 
be chosen to be mutually 'orthogonal' {QkQi)p — '5fc./(Qfc)^- Considering ob- 
servables with a vanishing equilibrium expectation (A) ^ — 0, strict positivity 
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of the right-hand side (RHS) of Mazur inequahty ([T|) is a convenient indicator 
of non-ergodicity of the observable A. Mazur has shown that in classical statis- 
tical physics, the inequality (jlj is merely a corollary of the Khinchin theorem 
[10) . Later, Suzuki jT7] has proven a quantum version of the bound ([T|), strictly 
applying only to finite quantum systems as his proof is based on explicit diagonal- 
ization of the Haniiltonian operator. Existence of non-trivial constants of motion 
which is characteristic of completely integrable systems implies non-ergodicity of 
(almost all/generic) observables, making an intimate connection between non- 
ergodicity and complete integrability, both in classical and quantum statistical 
mechanics. 

The inequality ([T]) has found numerous and very useful applications in con- 
densed matter physics as it is naturally suited for bounding dynamical sus- 
ceptibilities within the linear response theory. For example, within the frame- 
work of Kubo's linear response approach, the zero-frequency Drude peak (see 
e.g.[9 ) is defined in terms of the real part of (heat/electric/spin) conductivity 
(T^(cli) — 2'KDp5{uj) + a^^^{Lo). The constant Dp known as the Drude weight can 

be expressecF] for a one-dimensional quantum lattice of size n as 



The symmetrized correlation function is used in order to render Drude weight 
manifestly real. Here, in ([2]), J„ = 3^ is an extensive (energy /particle/spin) 

current operator, and is a current density at site x. Positivity of the Drude 
weight Dp > is a signature of a ballistic transport at finite temperature and is 
generically related to complete integrability via Mazur inequality ([T]) , as pointed 
out in Ref. 

It should be emphasized that the order of limits in the Kubo-type formula 
([2]) is crucial, namely taking them in a different order can sometimes produce 
completely different results (see e.g. Ref. [16J). Within Suzuki's approach, ther- 
modynamic limit (TL) of letting the system size n, or the number of degrees 
of freedom, n — > 00, can only be taken at the end, which is not according to 
a fundamental principle of statistical mechanics which states that TL n — >■ cx) 
has to be always considered before the long-time limit i — > 00 (say as in formula 
([2|). Our program in this paper is then to use the natural language of quan- 
tum statistical mechanics of systems with local interactions — the quasi-local 
C*-algebras — and to develop an algebraic approach to Drude weight bounds 
of the Mazur-type pertaining to strictly infinite quantum lattice systems. In the 
Suzuki's theorem conservation laws have to strictly commute with the Haniilto- 
nian, [H, Qk\ = 0, for all finite sizes n. For that reason, one has to often study 
periodic boundary conditions which may only guarantee existence of such exact 
conserved quantities. However, in our C*-algebraic setup we can accommodate 
also for quite common situations, where for any finite system size n the conserva- 
tion laws Qk might not be perfect, but the commutator dQfc/dt — i[H, Qk] may 
result in terms supported at the system's boundaries. In this way, TL n — >■ 00 
can be taken in the beginning, by inclusion of larger and larger open lattices, 
and no resort to periodic boundary conditions is needed. 

^ For a deta iled discussion and derivation of the linear response expression of Drude weight 
see subsection 16. 11 
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1.2. Summary of the main results. The main idea behind our susceptibihty 
bounds is to bring together two classical results in quantum spin lattice sys- 
tems with local interactions |5|[H1[T1 , 2 , namely (i) the exponential clustering 
property of finite temperature Gibbs states in one-dimensional translationally in- 
variant lattices, and (ii) Lieb-Robinson kinematic bound on the spatio-temporal 
propagation of local quantum correlations which result from boundary terms of 
the commutators. The projection form of the Lieb-Robinson bound [4 enables 
us to nicely separate the causality light-cone of a time-evolved observable in 
the autocorrelator and to exploit the property of exponential clustering in the 
equilibrium state. 

After introducing the notation and outlining the main concepts ((i) and (ii)) 
in section[2]we write and prove in section|3]our main result (Theorem 1): Namely, 
given any observable J„ and Hamiltonian of a finite lattice system of size n, 
such that the interactions are translationally invariant, and another quasi-local 
extensive observable Qn which has a property that the commutator Qn] — 
hi — 6„ where bi and 6„ are two operators supported at the left and the right 
boundary of the chain, we show that the Drude weight ([2| is strictly bounded by 
asimple expression Djj > ^ lim„_j.oo "■"'^(^(•AiQn + QnJn) pY / {Q'ii) j^- In section 
|4] we then provide trivial generalization of the result (Theorem 2) to the case 
where we have an arbitrary set of almost-conserved quasi-local operators. In 
section |5] we describe a nontrivial application of our results for bounding the 
spin Drude weight in the anisotropic Heisenberg XX Z model. In section |6] we 
conclude by discussing the assumptions needed to equate the thermal-averaged 
correlator ^ with the canonical Kubo-Mori expression, and state some remarks 
on possible other general contexts where results of our type may appear. 



2. Preliminaries 

We consider the following setup, where the notation essentially follows Bratteli 
and Robinson P]. Let N S Z+ be a local Hilbert space dimension, say N — 2 
for spins 1/2 or qubits, and 'Ql^ = C^^^ , x € Z, a local on-site matrix algebra. 
We associate a matrix algebra to any finite open lattice of integers, the so-called 
chain [x, y] = {x, x + 1, . . . ,y — 1, y}, as ^[x.y] = ^^z=x and define the quasi- 
local (UHF) C*-algebra 2t = 2tz in terms of a closure of the limit by inclusion 
[x,y] Z. 

Let the interaction h G 2t[o,d^-i] be an element of a local spin C*-algebra on 
dh sites 1^ The Hamiltonian 

n—dh+l 
x=l 

is an operator acting on a finite chain yl„ = [l,n], which is a sum of local 
energy densities = rj^ih) G Slj^j^i^-i-d^-i] , obtained by a group of lattice {shift) 
^-automorphisms rjx of 2t, defined by rjy{ax) = a.T+y Translationally invariant 



^ We shall assume > 2, as the on-site interaction with = 1 represent a trivial case 
with strictly local dynamics. 
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Hamiltonian can be understood in terms of the limit by inclusion yl — > Z of 
Hamiltonians Hj\ for arbitrary chains A 

max A~dh + 1 

Ha^ J2 (4) 

a:— min A 

The latter defines another group of ^-automorphisms of the quasi-local algebra 
21, namely the time automorphism 

n{a) = ]im T,-^(a), r^[a) := e'^-'ae"'^-*, (5) 

and a finite temperature equilibrium expectation, namely the infinite volume 
Gibbs state 

ujg{a) = lim — -. — , (o) 
' A-^z tr(e-'^^'i) ^ ' 

strictly defined only for local operators a and extended to 21 by continuity. Araki 
(Theorem 2.3 of [2]) has shown that such Gibbs state is an extremal (t, /3)-KMS 
state, which is invariant under space and time translations 

w^(?7x(A)) = ujg{A), (7) 
= ^^(A), (8) 

for any A S 21, x S Z, t G M. Most importantly, Gibbs state has an exponential 
clustering property (ECP) (Theorem 2.3 of [2 , see also Theorem 3 of [12J): for 
any pair of local operators / G ^[-df,-i], 9 G ^[o.dg-i], df,dg G Z+, and a 
displacement x G 2+ one has 

< ^ll/llllslle-"^ (9) 

where k, p are two positive constants, which do not depend on x, neither on /, g. 
Next, we define an extensive current operator in the open chain 

n—dj+l 
x=l 

where the current density j belongs to a dj-site local algebra j G 2t[o,dj-i]- Iii 
fact, the observable J may not necessarily be interpreted as a physical current, 
but it can be any spatial sum of a local self-adjoint operator j (representing an 
extensive translationally invariant observable) , the only condition being that its 
local equilibrium expectation vanishes 

^p{j)=0. (11) 

We shall think of it as a current merely because the most important applica- 
tion we have in mind is in the quantum transport. 



Definition 1. The key concept in our work is a quasi-local translationally in- 
variant conservation law Q with the following properties: 
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1. Q is a translationally invariant spatial sum of exponentially localized (quasi- 
local) operators. For any finite chain yl„; 



n-d+l 



Qa^^Y^QaI QaI^ ^i'^' <7exp(-Cd), (12) 



d=l 



where q'^'^^ — {q^'^'^)* <E ^[o.d-i], O'f^d o,^^ positive n-independent constants. 

2. Q should have vanishing thermal expectation value 

MQaJ^O. (13) 
We can thus assume also that all orders of local density g'^'^-' satisfy 

Mq^"^) = 0. (14) 

3. The operator Qa„ is almost-conserved on any open chain A„, i.e. it commutes 
with the Hamiltonian Ha,^ except for terms that are supported at the boundary 
of the chain 

[Ha^,QaJ^B9^ (15) 
where d„ = [1, dh] U [n — dt, + 1, n] and 

Bd„ bi - bn-d^+i (16) 

for some local operator h £ ^[o^di-i]- 

For concrete, nontrivial examples of Q, see section [5] 

Take now any chain A which is sufficiently bigger than yl„, say yl 13 [—dh + 
2, n + d;i — 1]. Then, the almost-coniniutation identity (|15[) impHes 



[Ha, QaJ - Ba„ + [h^, QaJ + [hn, QaJ (17) 

where 

n 
x— — dhA-2 x—n—dh-\-2 



represent the left and the right near-boundary interactions. The RHS of ( 17l can 
be rewritten as a sum of two quasi-local operators localized near the boundary 

[Ha,QaJ ^ Bi^ + Bn, (19) 



71 

Bn:=J2bK\ (20) 

d=Q 

x+dh-1 



n x — d-\-l 

^(0) ... hid) 



6f b'i':= E E [^-4'^]- (21) 

x—n—dh-\-'2 y—n—d—dh^^ 
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Note that the supportf|^of the boundary operators are 
supp b^^^ C [—dh + 2, max{(i;i — 2 + d, db}] , 

supp fcj^'' C [min{n — o? — d/j + 3, n — c?b + 1}, n + — 1], (22) 

and that they are exponentially localized (following from the definition (12 1 and 
elementary inequalities, ||A-B|| < ||v4||||i3|| and the triangular inequality): 

II^L'^II - < max{4(d. - m\hUme-^'- (23) 



Eq. ( 19 ) results in the Heisenberg equation of motion for the almost-conserved 
operator 

(d/di)T,^(Q^J = \[Ha.t^{Qa^)] = ir/i(i?L + i?R), (24) 

which together with the initial condition Tt=o(Qyl„) — Qa^-, S'fter taking the 
limit yl — >■ Z, integrates to an explicit time dependence 



n(OyiJ = Qyi„+i / dsr,(BL + BR). (25) 

Another crucial technical tool that we shall facilitate is the Lieh-Rohinson 
estimate (LRE) |11| which bounds the speed at which a disturbance propagates 
through a quantum spin system with local interactions. Let / g 2lx, g G 21;-, 
where X, _r C Z are two subsets with \X\, \r\ sites, such that at least one of 
them is finite. A useful form of LRE (see e.g. |3]) then states 

||[n(/),5]|| <</>min{|X|,|r|}||/||||5||cxp(-/i(dist(X,r)-«|t|)), (26) 

where dist(X, r') = rain^izx.yi^r \x ~ y\ is the distance between the sets X,r, 
and (f), fi and v are some positive constants, independent of /, g, and t. 

Regarding any finite subset F d we define a projection mapping (•)r '■ 
2t ^ 2t as 

{A)r ]im^ ^ ^^^^ " ITz/ MUA\r)UA\rAUA^r (27) 

where trx denotes a partial trace with respect to a local algebra supported on 



X. The RHS of (27 1 provides a very useful identity, where dii{Ux) denotes the 
normalized Haar measure for the integration over the full unitary group of a 
A^l^l dimensional Hilbert space over sites X. 



Bravyi, Hastings and Verstraete (BHV) [4 have reformulated LRE (26) in a 



very convenient way, namely they have shown that ( |26p implies the estimate 

lln(/)-(Tt(/))r|| < 01^1 11/11 cxp(-Mdist(X,Z\r)-«|t|)) (28) 



^ Due to existence of the principle of locality in the systems with finite-range interactions, 
the notion of an operator support enters naturally into the discussion. The support supp A of 
an observale A £ 21^ is the minimal set F G A for which A = A(X) lA\r foi" some A G 21 j^. 
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where / e 2lx, and r' C Z an arbitrary set of sites. The constants 0, /x, v are the 
same as in (26 1 and hence also do not depend on F. The BHV inequahty (28 1 



is proven by applying the identity ( 27 ) , writing 



Wnif) - (Tt(/))r|| = il ]im / d^^{U^\^)[rtif),UA\^]UX\^\ 



<\irn I dfi{UA\rmUA\r.Ttif)]l 



(29) 



and then using LRE (|26|). Note that C/*^^ = U^^ 



A\r- 



3. Thermodynamic limit of Mazur inequality 



Let us consider a suitable extensive translationally invariant observable J (10 1 



Then we define the Drude weight in terms of a time-averaged autocorrelation 
function as the following double limit, 



Dp :— lim lim 



2n 2t 



dt'uJi3{JA„Tt'{JAj)- 



(30) 



In order to avoid any ambiguity and to make the definition precise, we stress 
that the Gibbs state utp and the time-evolution Tt on the RHS of (30) are already 
taken for an infinite system, before the infinite volume limit n — > c» is applied 
to the extensive observable Jj\^ . However, we will show later that RHS can be 
expressed in terms of local quantities only, therefore a single TL suffices. 

The Drude weight Dp is an important quantity in linear response theory of 
condensed matter physics (see the discussion in subsection 6.1 1. In the following 
we prove a general and useful inequality related to it, and clarify its existence 
in infinite one-dimensional systems. Moreover, Dp can be considered as an in- 
teresting ergodicity indicator of C* dynamical systems. 



Theorem 1. (i) The following limit exists for any t G 

n 

c(t) = lim V ujp{j^Tt{j)). 



(31) 



(ii) Let us assume that the symmetric time-average of c(t) exists 

1 



lim 

i— voo 2i 



ds c(s). 



(32) 



together with a very mild condition, namely that the integral of \t\{c{t) — c) grows 
slower than 



lim 

t-i-oo t^ 



ds\s\{c{s) - c) = 0. 



Then, the double-limit Dp {30) exists and is equal to 



(33) 



(34) 
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(Hi) Let Qa„ be a self-adjoint almost-conserved quantity, satisfying Definition IT] 
(Eq. 15). Then, under the assumptions of (ii), a lower hound on D/j exists ana 
is equal to the limit 



2 Ti-i-oo n 4w/3(Q^_J 



D0>^ lim 3 ' . (35) 



{A, B} = AB + BA denotes the anti-commutator. 

Proof. We start by considering the following finite-time-averaged self-adjoint 
operator 

AA„,t-=^l f dt'{Tt,{JA,J-aQA,J, (36) 



y/n t 

where a G M is a free parameter. Since the state up is a positive linear functional, 
we have ujj3{A\ t) > for any t,Q; G M, n G Z+, or equivalently: 

0<^ f dt' f dt"lLOp{Tt>iJAjTt„iJAj) (37) 



Jo 



n 



-jf dt'-{0Jfi{Tt'{JAjQAj+MQA^Tt'{JAj)} (38) 

t Jq n 

+ —upiQlJ. (39) 
n 

We shall proceed to show that TL n — ^ oo of all the terms in the above inequality 
exists, treating it term by term. 

Let us first discuss in detail the mixed terms ( 38 1 . Time invariance of the 
Gibbs state Q implies that the integrand of (38 1 can be rewritten as (writing 
now time integration variable as t) 

UJf3{JA„T^t{QAj) + UJ(i{T^t{QAjJAj- (40) 

As both terms can be treated on equal footing, we shall focus on the first one 
u)p{JAnT^t{QA„)) a-nd show that in the limit n — > oo the term becomes time- 
independent. Using the explicit time evolution (25 1, the linearity of the state w^, 
and an elementary integral triangular inequality, we estimate 

|W/3(JA„T_f((9yl„)) - Ui3{Ja^QaJ\ 



< 



J ds \u;p{Ja„Ts{Bi^ + Br))\ 



-t 

<l ds {\ujpiJA„T,{Bi^))\ + \u;p{JA^T,iBnm. (41) 



We note that, before taking the TL n — > oo, the mixed terms (38 1 have to be 
multiplied by 1/n, hence it is enough to show that all terms on the RHS of (41 1, 
for arbitrary time variables s, do not grow with n, i.e. are bounded by constants 
that do not depend on n. Again, it is enough to focus just on one of the terms in 
the integrand, say \u}i3{ Ja„Ts{Bi,))\, while the other one can be treated in exactly 
the same way. Expressing the operators in terms of local densities, we obtain: 
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Fig. 1. A scheme visualizing the estimation of the crucial term |cj^{ja;Ts(6]^ ))|. The operator 

Tsih^"'), supported near the left boundary at time s = 0, is localized within an effective 
light-cone (light-green region) and can be approximated with an error exponentially small in 
half-gap i (indicated) to the support oi j^. (red strip), by projecting it onto a sublattice F 
(painted in yellow). 



|c.MJ^„r,(i?L))l < E E (42) 

For the time being, let us fix the integration time variable s, and the summation 
variables d, x. We first assume that the Lieb-Robinson light-cone emitting from 
the local operator , at time s. is not including x, i.e. there is a positive 
half-gap £ > 0, defined as (see Fig. [T]): 

l^\^{x~ v\s\ - max{o?,j - 2 + d, 4}) . (43) 

Defining the chain of sites including the light-cone plus the half-gap of sites on 
each side as 

r=\-dh^2- [w|s|-f-£J,max{4-2-f b|s|+^J] (44) 



we can estimate the term on the RHS of ( 42 1 as 

\u^f,{j,,T,{hf ))\ < \u:piMMbi'^))r)\ + IMMrsibi'^) ~ (r.(6L'^))r))|, (45) 
where the first term is further estimated using the ECP (|9|, also noting (111, 

|c^Mj.(T.(&[''))r)| - \u^,[3.[r,{^t^))r)-^ii{3.)^iAkrsi^t^))r)\ 

<A\3\\\\{rs{h'f))r\\er^'^'-^\ (46) 



while the second term is bounded by BHV form of the LRE ( 28 1 

\^,{3.{rs{h'^l^) - {rs[ht^))r))\ < \\j\\\\rsibi'^) ~ (r.(6[''^))r|| (47) 
< 0max{dfc + - 1, 2d,, - 3 + d} max{d,,(d,, - l)-f\\h\\,\\b\\}\\j\\e-^'^'-^'>-^'' . 

The norm of projected evolution can be estimated generously by another use of 
the triangular and the BHV inequalities 

UrAbi'^))r\\ < \\rsibi'^)\\ + \\rAb['^) - {rAbi'^))r\\ (48) 
< (1 + 0max{4 + dh-l, 2dh - 3 + d}) 
xmax{dft(d,,- l)7||;i||,||6||}e-«'*. 
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Putting Eqs. ( 46|47|48 1 together, writing the minimal exponent A = minjp, /i} > 
0; and taking suitable constants C, C > 0, which only depend on the local 
operator norms and dimensions but not on the size n neither on the variables x, 
d, and s, we have: 



(49) 



For the bound to remain valid for an overlapping light-cone, < 0, we should 
simply make sure that we choose C and C" large enough to satisfy the naive 
bound \ujp{jxTs(b^i^^))\ < |lj|l||6L'^'*|l. The whole term (42) can now be estimated 
as, introducing k{d) :— lv\s\ + max{dh — 2 + d, db}\: 



n—dj+l 



oo 

< "^{Cd + C')e-^'^ 



inax{0,n— 



k{d) 



E 



x'=0 

db^2 



d=0 



E' 

x'=0 



-\x'/2 



Y,{Cd^ + C"d + C"')e-^'' 
e«((e« + l)C + (e«-l)(C" + (e« 



1)C"")) 



(e« - 1)3 



=: K < oo. 



(50) 



where C" = C + C{v\s\+dh+db-2 + {l-e-^/'^)-^) and C" = C {v\s\+dh+db- 
2+ (1 — e^"^/^)^^). As C" and C" are at most linear functions of time s, we have 
also that K = K' + K"\s\ where K' and K" are n and s independent constants. 
Exactly the same estimates applies for the other term \uji3{Ja„Ts{Bji))\, so we 
have finally shown that the difference on the left-hand-side (LHS) of (41 1 is 
bounded by 



luJfiiJA^T^tiQAj) - C^M^A„Q^JI < ^K'\t\ + K"t\ 

We shall now study the convergence properties of the sequences 

1 



(51) 



(52) 



and show that they are, in fact, Cauchy sequences. 

Let us consider, for the time being, an abstract sequence of this type 



1 



(53) 
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where 

n—df + l n — dg+l 



with ujp{f) = (^p{g) — 0. Exploiting the translational invariance of the Gibbs 
state ujpifxgy) = ujpifgy^x) we find (see Fig. [2]) 



n—df + l n—dn + l 



n 

x—l y— 1 

max{0.d/— dg} 
r—min{0, df —dg} 

n — dg 

+ E (i-^)-,(/ff.) 

r— max{ 1 — (ig + 1 } 
inin{ — l.d/— dg — 1} 

+ E (55) 

r— — n+df 

Clearly, ECP (|9]), combined with the trivial norm-bound, 

\i^{fgr)\ < ||/||||.9||min{l,KCxp(-p(r-(i/ + l)),Kexp(-p(-r-dg + l))} (56) 
guarantees boundedness of the sequence \v„\ < V < 00. But we have more: 



\Vn+l-Vn\ 



n{n + 1) 



max{0,d/— dg} 

(max{d/,dg} - 1) E ^pUgr) 

r— min{0,c/ / — c/g } 
n+1 — min{ — dg — 1} 

+ E (c^s - 1 + ''V/3(/5r) + E {df r)u}p{fgr) 

r=max{l,d/— + r= — ri — 

II/IIII5I 



< 



-{(max{d/,dj-l)(|d/-d<,| + l) 



n(n H- 1) 

df—l n+1— rfg 

+ E (rfg-l + + « E (rfg - 1 + Oe"''^''"''^^'^ 

r— max{l,(i/— dg + 1} 

miii{ — (ig — 1} ~f^g 

+ E ('^Z - 1 - ^) + E (^/ - 1 ~ r-)e-''(-'-'^.'+i) } (57) 

r— — dg-\-l r— — n—l-{-df 

WfWM 

n{n+ 1) 



|4 + c^g + rf/c^s - 2d/ - 2dg + 1 + 2k E (^+^/+^s-2)e"'"'| 



where we have adopted a convention that X]r=a;(---) = if x > Finally, 
we complete the sum in the curly bracket to a geometric series and simplify the 
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y 




1 ^^ii~^f '^f! n-df + l 



Fig . 2 . A scheme visualizing the double summation over x and y for computing the sequence 
Vn \55\ . Yellow, green, an d re d region, represent the terms under the first, the second, and the 
third sum on the RHS of | |55[ , respectively. Light-green and light-red stripes denote the te rms 
which are separated out as tne first (left) summations in the last two lines of expression l |57| l 
before ECP (loj is applied to the remaining terms. 



denominator to arrive at the quickly decreasing bound on the difference between 
the adjacent terms 

\Vn+l - Vn\ < +4+ djdg + V{df + - 1) + C}, (58) 

where v :— 2{ ^J^_-^ — 1), ^ '■— (ep^i)'' ~ ^- "^tiich proves that {w„} is a Cauchy 
sequence, i.e. hm„_>.oo exists and is finite. In fact it is equal to the limit 



lim Vn = lim ivpifg^;), 



(59) 



since it can be shown - again using ( 55156 1 - that the difference of the terms of 
the sequences on the LHS and the RHS of (59 1 is bounded by const/n. 



Similarly we now treat the sequence w„ = n ^a;^( Jyi„Q/i„)- Writing Wn = 
n~^'^^^iUJl3{JA„Q^A^), and associating f = j , g = q^'^\ we arrive at (again in 
the last step completing the geometric series) 



\Wn+l -Wn\<^J2 ^"*''(^' + ('^J + + 4 + "("^1 - 1) + C) 



< 



7||j|| 2 + (e^ - l)(d, + 3 + ^) + (e^ - l)^((d,- + 1 + i^)d, + 1 + Q 
n2 (e« - 1)3 



(60) 



and with the same bound for the 'transposed' sequence = n ^ujp{Qj\^Jj\^). 
Analogously, for the sequence u„ = n^^ujp{QA„QA„) = J2d,d'=i ^I3{Q'^aIQ^a]) 
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we obtain 



+i-M„|<\ y e-^^''+'''Hd^ + d'^ + dd' + i^{d + d'-l) + 



'd' = l 



< 



72 C - + e«(2 - 2C + e«(C + + 3)) 



(e«-l)4 



(61) 



We have thus shown that the hniits 



w = hni w„ = V V ujpijqi'^^), 

00 

n— ^00 -"^^ — ^ -'^^ — ^ 

00 

u= lim u^= 5]a;^(gWgf)), 



(62) 



i,d'=i xei 



exist and are finite (see Eq. (59 1 for the RHS). Moreover, Eq. (51 1 impUes that 

(63) 



the distance between w„ + w„ and the integrand in ( 38 1 , 

Zn{t) := n^^ {ujpiTtiJAjQAj +(^l3{QA„Tt{JAj)} 

is decreasing in n 

4:K'\t\ + 2K"t^ 



\Wn +w'„~' Zn{t)\ < 



(64) 



so the two sequences should have the same limit, i.e. lim„_j.oo Zn{t) — w + w' . 
Therefore, the TL of (38 1 exists, it is independent of t, and equals w + w' . 



At last, let us devote also some attention to the first term ( p37| ), which af- 
ter exploiting the time-invariance of the Gibbs state and substituting for the 
integration variable s = t" — t' , reads: 



1 



Unit) = - / ds 1 - — C„(s) 



(65) 



where 



Cn{t) 



n—dj 

^-MJaMJaJ)^ E (l-^^^)wM>n(j)). (66) 



r— — n-\-dj 



For each fixed time we can again show that {cn{t)} is a Cauchy sequence, 
namely 



^ 00 

K+l{t) - Cn{t)\ < ^^^^^^ J2 (1^1 + - ^)MrTt{m < 



L\t\+L' 



(67) 



r— — 00 



where the constants L, L' follow from bounding the sum of exponentially decay- 
ing (in a;) envelope of the spatio-temporal autocorrelation function 



MjMm < ||.7!pmin{l,e-^(l-l-'^^-''l*l)/2}. 



(68) 
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Eq. (681 can be proven in exactly the same way as Eq.(49l - combining the 
BHV inequahty and the ECP - simply taking j to substitute the boundary 
operator This implies that the limit of Cn{t) exists. It also coincides with 
the expression (31 1 of c{t), as the difference between the terms of sequences (66 1 
and ( [31] ) can be bounded by const /n (again using (68 1). 
This proves the point (i) of the theorem. 

As Cn(t) converges uniformly in t, the limit a{t) = lim„_i.oo an{t) of (65 1 exists 
as well and is equal to 



ait) 



1 



c{s). 



(69) 



Similarly, a uniform convergence of c„ (t) can be used to express the Drude weight 
([30| as 

rt a 



D 



13 



^ lim — lini 



dsc„(s) 



^ lim — 
2 t-s-oo 2t 



ds c(s) 



P. 



(70) 



proving the point (ii) of the theorem. 

Finally, writing the TL n — > cx) of all the terms in the inequality ( 37||39[ ) we 
have 

Q<a{t)-a{w + w')+a'^u. (71) 
One should observe an obvious identity 



a{t) - c = 2 



ds c(s) — c 



ds|s|(c(s) - c), 



(72) 



where, according to the assumptions ( |32[ ) and (33), the limits t 
terms on the RHS exist and vanish, hence lim4_j.oo a(0 



00 of both 
c. Taking then the 



limit i — >■ cx) of the inequality ( 71 ) and optimizing it with respect to the free 
parameter a, we arrive to the final bound 



c > 



w') 



l\2 



(73) 



This proves the point (iii) of the theorem. 

Let us conclude this section by making a few remarks. 

Remark 1 . The key in our proof was a convenient form of the ansatz ( 36 1 . It has 
been inspired by combining a finite-time-average version of the original idea of 
Mazur |r3] with an expected central limit theorem behavior in size n for local 
spin chains (hence the ^1 \/n prefactor) which allows taking the TL n 00 first. 

Remark 2. It is not clear at present if in the context of quasi-local C* dynamical 
systems, the assumptions of point (ii) of the Theorem 1, e.g. on the existence of 
the time average c, are in fact needed, or can be separately proven. In cases of 
more general dynamics it is of course easy to come up with counterexamples of 
dynamics for which c does not exist. For example, one may formally construct 
a gaussian process with the (2-point) correlator c{t) = ( — 1) Liog2(*/to)J fQj. -^Jiich, 
clearly, the average c does not exist. 



Remark 3. Note that the RHS of our estimate ( 35 1 can in fact vanish in some 
cases, e.g. when the limits w, w' vanish. In such cases our bound is nearly trivial, 
and merely expresses the fact that the Drude weight as defined by ( 30 1 is always 
non-negative. 
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4. Generalization to the case of several almost-conserved quantities 

The bound (|35| of the Theorem [l] can be readily generahzed to the case where 
one has several almost-conserved quantities. 

Theorem 2. Let (5[fc]A„; fc = 1,...,to 6e a set of m self-adjoint, almost-con- 
served quantities with densities q''^^^ G ^io,d-i]j ^,11 satisfying Definition jlj (Eq. 
I j^p . Then, under the assumptions of the Theorem 1, together with the points (i) 
and (ii) of Theorem 1 we have the following lower hound on Dp (30) 



iwi 



k,l=l 



where 



Wk 



1 1 °° 



(74) 



(75) 



= 1 x& 



and U ^ is an inverse of m x m covariance matrix 

1 ^ oo ^ 



d,d' — \ x£l 



Proof. The proof involves exactly the same steps as the proof of the Theorem [T] 
except that the ansatz ( 36 1 is replaced by a more general one 



1 1 

y/n t 



J dt' ^t'{JAj - ^a/cl9[fe]/l„^ , 



(77) 



with m real free parameters . The final step of optimization of m-dimensional 
quadratic form (in ak) then results in (74 1. T he lim it-identity (59 1 is used to 
write the compact expressions on the RHSs of (75|76l. 



5. Examples 

As a pool of nontrivial examples let us discuss the anisotropic Heisenberg spin 
1/2 chain (the so-called XXZ model). Here N = 2, and the local algebra 2l[o] is 
spanned by Pauli matrices cr*, s € {0,x,y,z}, with — 1, or cr* = ^{a^ ztia^). 
The Hamiltonian density reads h = aQdf + <Jq(jI + Aa-Qa^, dh — 2, where A is 
the anisotropy parameter. 

Grabowski and Mathieu [SI have shown that an infinite sequence of nontrivial 
local conservation laws of an infinite XXZ chain can be formally constructed 
using the boost operator B = ^ 'l2xe2^^^^ namely 

Q[k+i] = [B,Q[k]], k^2,3,..., Qi2]=Y.'^-' (^8) 
supplemented with the trivial on-site conservation law - the total magnetization 
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Translationally invariant operators (|78| which can clearly be written in terms 



of fc-site densities Q^^] = Y.x(^z'nx{q[k\h q[k] e 
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[0,fc-l] 



are strictly conserved, 



I.e. they exactly commute with the Hamiltonian Qp] among each other 
[Q[fc]jQ[;]] — 0, only for infinite lattice (Z) or for periodic boundary conditions. 
Note however, that they do not exist as elements of C*-algebra 2t in the former 
case, while in the latter case the 'limit by inclusion' construction yl — Z does 
not work. 



In the setup of C*-algebraic statistical mechanics, the operators (78) should 
be considered on a finite open chain yl„, written as Q[k]An^ being the elements 
of 2t. Interestingly, it has been shown again by Grabowski and Mathieu [5] 
that for open boundary conditions, i.e. considering the Hamiltonian Ha^, half 
of the conservation laws (for odd k) are destroyed, while the other half (for 
even k) can be amended by adding terms supported near the boundary of 



[2l]A, 



\2l' 



]A„ 



0, for 



An- Q'[2l]A„ = QmA„ + Ql2l]L + (9[2/]R, SUCh that [Qj, 

1,1' = 1, ... , [n/2j . Nevertheless, as we have shown in the sect _ 
need exact conservation to state dynamic susceptibility bounds( 35|74 1. There- 



fore, we rewrite the procedure (78 1 of 
and boundary remainders pj^j G !<!l[o.fe] 
relations 



in terms of densities 



one does not 



qik] 



e 21 



[0,k-l] 



such that one has almost- commutation 



'^[HA„,Q[k]Aj = Vl{P[k]) - Vn~k{P[k])- 



(79) 



Operators which can be interpreted as the current densities, are actually 
determined from the local continuity equations (d/dt)g[/c] = r]-i{p[k]) — P[k], or 



Plk] 



mipik]) ^^^[hx,vi{qik])]- 



(80) 



Local-algebraic version of the boost relation ( 78 1 gives the other recurrence re- 



lation which determines the charge density in the next order in terms of the 
charge and current densities of the previous order: 



Qlk+l] 



-.Plk] 



k-1 



'-Y,{^ + ^)[hx,q[k]], fc = 2,3. 



x=0 



(81) 



Clearly, p[i] — j — 2(crQcr^ — crQcrf ) is just the spin current, g^] = h, while the first 
few higher charges and densities can easily be obtained solving the recurrence 
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( 80|81 1 by means of some computer algebra]^ 



yzx 



(jxyyx ^xzxz ^ ^yOyO _ ^yxxy _|_ ^yxyx ^yzyz ^ ^zOzO ^zxzx _|_ ^zyzy 

2^zzOO>) _ 2cr'<^^"'^ - 20-^^'^^, 

3 _u lA/W/rXyzOO _ ^yxzOO j_ ^zxyOO ^zyxOO^ ^ RA2|'_^x0yz0 _^ ^xyxxz 



i^xyyyz _|_ ^yOxzO _ ^yxxxz ^yxyyz _ ^zxOyO _^ ^zxxxy _ ^zxxyx ^zxzyz 
_,_ ^zyOxO _^ ^zyyxy _ ^zyyyx _^ ^zyzxz-, ^ ^^q^2 _^ 8)(_^xzyOO _^ ^yzxOO) 
+ QA{a^^^^'^ — a^y'^^'^ — ^-xy^zx _ ^xyyzy _|_ ^xzOyO _ ^xzxxy _|_ ^xzxyx _|_ ^xzzyz 
_ ^yOzxO _|_ ^yxOzO _|_ ^yxxzx _|_ ^yxyzy _ ^yzOxO _ ^yzyxy _j_ ^yzyyx _ ^yzzxz 
_ ^zOxyO _|_ ^zOyxO ^zxzzy _ ^zyzzx-j _ g^xzzzy g^yzzzx 

P[2] = -29[3], 

P[3] 2A^{<7^''''^ + o-^yy") - (2z\2 + 2) ((7°'^'^" + cr°yy") + 2z\(-2(7"""" - d'^y'^y 

_|_ g-xyyx _ ^xzxz ^yxxy _ ^yxyx _ ^yzyz _ ^zxzx _ ^zyzy^ _|_ 2a^™^ + 2a^™^ 
P[4] = (4Z\3 + 4Z\)(cr°''y^° - Cr°y''^° - fJ^'^y™ + a^y'^^O) + 4Z\2(_crXyxxz _ ^xyyyz 



_|_ ^yxxxz ^yxyyz _|_ ^zxOyO _ ^zxxxy _|_ ^zxxyx ^_ ^zxzyz _ ^zyOxO _ ^zyyxy 
_|_ g-zyyyx _ ^^zyzxz-j _^ (4/\2 _|_ 4)( — o-"''^^^ + cr^-^^'^^ + (jXzyOO _ ^yzxOO-j 
+ 4Z\(2(T°™y" — 20'"^^'^° + u^y'^^^ + 0-xyxzx _|_ ^xyyzy _ 2fjXyzOO _ ^xzOyO 
^xzxxy _ ^xzxyx _ ^xzzyz _ ^yxOzO _ ^yxxzx _ ^yxyzy _j_ 2jjyxzOO _|_ ^yzOxO 
g-yzyxy _ ^yzyyx _|_ ^yzzxz _ ^zxzzy _|_ ^zyzzx-j _|_ ^^xzzzy _ ^^^yzzzx ^g2) 

where we use the notation cr*i'*2- -«fc = cr^i (g) cr'^a ^ . . .(j^k^ g {0,x, y, z, ±}. It 
has been shown in Ref. |18j that Q^^ can be applied to the Mazur inequality 
|131I17| with periodic boundary conditions in mind, to yield finite, non- vanishing 
spin Drude weights for a general XXZ model with a transverse magnetic field 
of strength x added to the Hamiltonian densitjj^ h' = h + xcg, at finite (non- 
zero) or even infinite temperature (where the evaluation of expectations uj 13(A) 
becomes easiest). This interesting result called for deeper theoretical understand- 
ing, as the Suzuki's proof [17_ only allows to consider the incorrect order of limits, 
namely i — > 00 first, for a finite periodic system, and only then n — > 00. This 
problem has been now settled in the present paper. 

Furthermore, it has been pointed out in [18J, that for a vanishing external 
field X = 0, the Mazu r bound on the spin Drude weight always vanishes, for any 



set of from (78 1. This results from opposite 'spin-fiip' symmetry 5 : tr^ — >■ 
— ct^jCT''' — > cr^ of the spin-current, Sj — —jS, and of the conserved charges, 
^Qlk] = Q[k]S. Non- vanishing magnetic field breaks the spin-flip symmetry of the 



Hamiltonian density and makes it possible that coefficients Wk ( 75 1 are non- 



vanishing and the bound (74) is strictly positive. On the other hand, there has 



been clear numerical evidence (see e.g. |S]) suggesting finite Drude weight and 



* Using the Mathematica code http : //chaos . f mf .uni-lj . si/prosen?action=AttachFile&do 
=get&target=PauliAlgebra.nb one may obtain explicit form of gjj.] and P[k-i] up to fc = 10. 

^ Note that considering non-zero transverse field x is - in TL - equivalent to considering 
symmetry sectors with non-vanishing total (conserved) magnetization. 
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ballistic spin transport at any temperature even at vanishing magnetic field 
strength x = 0, in the easy-plane regime \A\ < 1. This seemed to suggest that 
another, nontrivial quasi-local conservation law should exist with the same spin- 
flip symmetry as that of a spin current. 

Indeed such a missing translationally invariant quasi-local conservation law 
Q with negative spin-flip symmetry has recently been found [TS], in case when 
|Z\| < 1, satisfying almost commutation condition 



[HA„,QA„] = -2ial + 2ial 



(83) 



It has been shown that Q admits a simple matrix-product representation in 
terms of infinite rank, almost-diagonal matrix operators Aq, A±. acting on an 
auxiliary Hilbert space with orthonormal basis labeled as {|L) , |R) , |1) , |2) , . . .}, 



An = |L) (L| 



|L)(1| 



A_ 



ID (Rl 



|R) (R|+^cos (r(^)|r) (r| , 

OO 
OO 



r=l 
2 



r+1 




2 





r) (r+l\ 



(84) 



+ l) ^) |r+l) (r|, 



where A = cosip. Namely, Q satisfies the conditions of the Definition 1, with 
local densities of order d> 2, (q*^^^ — 0), generated as 



L| A+A,, 



■ -A-s^.i A_ 



,(^) = i 

S2,...,sd_ie{0,±} 
For example, the first few orders read explicitly: 



|R) 



(85) 



r(-S2).--(-Sd-i) 



q(2) 


= i{a+- -cr- 


-+) = U 


q(3) 


= iZ\(a+"- - 


--"+), 




= iZ\2(^+oo- 


- a-""+) 



(86) 



-++ 



For the resonant values of the anisotropy A = cos(7r//m), for coprime l,m CzIj, 
TO > 1, the rank of matrices (84 1 becomes finite, i.e. it is m-f 1, whence the norm 
of q^'^^ can be estimated < je^^"^ by powering the following rank-(TO -I- 1) 

transfer matrix 



|L) (L| 

1 . 2 

- sm 









+ 1 


R) (R| 












r + 1 






2 


TO J 



|L)(1| 



1 R 



E 



COS 



\r) (r + l| + ^sin2[(2 [^J 



7rr/\ 

TO J 
TO 



(87) 



|r + l) (r| 



whose subleading eigenvalue is strictly smaller than 1 |I5| . In this regime, one can 
use (85 1 in our Theorem 1 to bound Dp even for a vanishing transverse magnetic 
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field X = 0. The calculation can be made exact for infinite temperature /3 — >■ 0, 
where we present a rigorous lower bound expressed as |15| 

Db 1 " 

lim > ADz, Dz := - lim 



/3^0 /? - " 4n^oo (L|T»|R) 

This expression can be evaluated explicitly [Ij in terms of Jordan decomposition 
of T, yielding a fractal (nowhere continuous) dependence on A: 

Dz = -{l-A^) Z\ = cos-. 89 

2 TO — 1 TO 

Note that two marginal cases with \A\ — 1, (for / = and I = m), are 
exceptional, as translationally invariant operator 

oo 

Q{A = ±1) = i^Z\'^-2^7y,(a+^- - a^aj) (90) 

is no longer a spatial sum of exponentially localized operators (therefore incom- 
patible with the Definition [I]) , i.e. ||(?''''-'|| ceases to decay with increasing d. 

Computing the Drude weight bounds for a finite temperature (/3 > 0) is 
certainly more tedious and cannot be done as explicitly as for the infinite tem- 



perature (88 1. One possible approach for high-temperature is, for example, a 
/3-expansion |7j . 

We stress that the bounds provided by our Theorems 1 and 2 for the XX Z 
model are mathematically rigorous if one assumes in addition the existence of 
the Drude weight. The second assumption in the Theorem 1 (ii), essentially re- 
quiring an (arbitrarily slow) relaxation of the extended current-current temporal 
correlations c{t) — c, can be understood as necessary for having a uniquely defined 



Drude weight (see the discussion in subsect. 6.11. 



6. Discussion 

6.1. On the linear response derivation of the Drude weight. Note that a slightly 



different form of the Drude weight than ( 30 1 follows from a strict derivation of the 

linear response where the canonical Kubo-Mori inner product dXuj[i{a*Ti\{b)) 
replaces a simple thermal average ujis{a*b). We will show that a small extra 
assumption of "non-ergodic dynamical mixing" is needed in order to justify the 
simple thermal-averaged expression ( |30p . Let us here carefully outline our linear 
response setup and the main steps of our argument. 

We consider the constant-gradient field perturbation to the Hamiltonian, 
which extends on a finite symmetric sublattice [— n,n] (for A D [— n, n]), as 

n — dg + l 

H^'''=Ha-F ^ xq.j;, for some q € 2l[o^d,-i]- (91) 

^A-i-i, generates a perturbed time evolution T^'^\a) = limyi_j.z e'^-^ ae~^^^ , 
on an infinite lattice, but still for a finite field extension [— Starting in the 
equilibrium state LOp for = 0, at i = 0, and then quenching the Hamiltonian 
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by switching on the force field (91 1, we define the canonical Drude weight as the 
asymptotic rate at which the local current in the bulk increases per unit time, 
after taking the limit of infinite field extension: 



Dr- 



um — lim 



dF 



(92) 



F=0 



Note that, in contrast to Ref. [5], we have to take a vanishing force limit F 
first, in order to make our perturbation bounded (and well defined) in the infinite 
extension limit n — >■ c». Writing in the first order of Born/Dyson expansion 



){3)=J-^Ff ds a;[r_,(g,),j]+0(F2), 



(93) 



and observing the following identity, for any f,gE'Ql 



dXu0{n.Mf))g) 



(94) 



where 5{f) = lim/i^z i[-ffA, /] is the ^-derivation with respect to the the unper- 
turbed dynamics, we arrive at 



Dr — lim — lim ds dA > x uja ( 



-iA iH<lx))j) 



(95) 



Finally, we choose the current j and the force-field density q, such that they 
satisfy the continuity equation j^-i — jx = S{qx) (see, as an example, Eq. (80 1). 
The sum over x can now be expressed as a difference of two sums which result 
in, after shifting the index x x + 1 in the first sum containing jx-i, 



1 n fp ^ 

Ds — lim — lim ds dA > ( 

P i^oo 2t «^oo J„ Jq ^ P"- 



-iX {jx)j) , 



(96) 



except for the boundary terms of magnitude \i^i3{T-s-ixUx)j)\ at a; w ±n which 
decay exponentially in n and uniformly in z = — s— iA (and thus do not contribute 
in the limit n — oo of ( [96| ) according to the Theorem 4.2 of Araki [2_: Namely, 
lim„_>ooel"l''||[/, T2(?7„(.g))]|| = 0, applying for any strictly local f,g and z £ 
C. Shifting the time evolution to the second factor of RHS (96), the sequence 
Cn{z) := (jxTzU)) couvcrgcs uniformly to c{z) (expression (31) of 

Theorem 1, but for complex time argument z), again as a consequence of the 
above mentioned Araki's proof of locality of complex time evolution. Thus, the 
limit n — > oo and s, A integrations can be interchanged. Finally, we obtain 



D, 



lim — 



ds 



d\c{s + iA). 



(97) 



Clearly, a sufficient condition for the equality between the thermal ( 34 1 and the 
canonical (97) Drude weights 

Dp - bp (98) 
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is a kind of non-ergodic weak-mixing condition 



1 /•* 

lim — / ds|c(s + iA) - c| = 0, VAg[0,/3]. (99) 
t->c» 2t Jq 



This follows from writing — Dp{t)\ < ^ ds dX\c{s + iA) — c| where 



Di}{t) and Df}{t) denote the expressions (30 1 and d97|, respectively, before taking 



the limit t — > oo. For example, for the condition ( |99| to be fulfilled it is enough 
that the correlation function c(z) converges (to an asymptotic value c, no mat- 
ter how slowly), as \z\ — > oo anywhere within the thermal strip Imz G [0,/3]. 
This is likely to be expected as a generic property of many-body systems with 
thermodynamically dense spectra, however it seems one may not trivialize it 
using Cauchy's integral theorem, as in a related study of Ref. [S] with holomor- 
phic dynamics of finite-system's observables. For instance, our c(z) may not be 
holomorphic (due to an infinite extension limit n — > oo needed for a definition 



(31 1). In fact, a better understanding of analytic properties of c(z) seems to be 



an interesting problem for future investigations. 



Nevertheless, even if Dp (30 1 may not be both, rigorously and generally con- 
nected to the linear response theory of the Drude weight, it might still be inter- 
esting general-purpose quantity in its own right, as a simple dynamical indicator 
of non-ergodic C* dynamical systems possessing non-trivial conservation laws. 



6.2. General remarks. Note that our Theorems Tj2 can be applied to any com- 
pletely integrable quantum chain, like the one-dimensional fermionic Hubbard 
model, supersymmetric t-J model, etc., for which one has algebraic procedures 
(e.g. algebraic Bethe Ansatz) for obtaining nontrivial conserved quantities which 
can be written as exponentially convergent sums of local operators. It is enough 
that these operators are conserved only in the TL, while for finite chains their 
time-derivative is supported only on finite domains near the boundaries of the 
chain. 

We conclude by remarking that the conditions of our theorems could be re- 
laxed or generalized in some cases: (i) The translational invariance of interactions 
was required only to guarantee a sufficiently strong clustering in the equilibrium 
state. This can perhaps be relaxed, and replaced by an algebraic clustering in 
some instances (with inverse power larger than 1, e.g. for zero temperature equi- 
librium states in the so-called critical systems) , or exponential clustering may set 
in for other reasons (e.g. due to disordered interactions and localization, however 
in such cases we are not aware of non-trivial examples of extensive conserved 
quantities), (ii) Quasi-local almost-conserved quantities need not be exponen- 
tially localized, but it would be enough to have a power law bound Hg*-''^ || < ^d~'^ 
for sufficiently large inverse power v > 3. (iii) Similar theory could be built for 
higher, ZJ-dimensional quantum spin lattices, where time-derivatives (commu- 
tators) of almost conserved quantities could result in terms distributed on the 
(D — l)-dimensional boundary of the lattice. Again, we do not see yet any ap- 
plication of this potentially very interesting generalization. 
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